400	441	80	1	! Ne, Nn, Nb, Nm
4	  2			! Nd, Ng
1.				! Sx
0	0	-1		! ipar(2:3), ipar(6)
1.e-9	1.e-9	0.		! fpar(1:2), fpar(11)
--------------------------------------------------------------------------------
! Conduction in an orthotropic rectangle												
! from	H.S. Carlslaw and J.C. Jaeger
! Conduction of Heat in Solids, 2nd. ed., Oxford, 1959, pp. 166-169											
! case V	(Eq. 26) p. 169											
! with specific choice: f(x) = sin (m pi x/a) so that	Am=1 for m=n and Am=0 for m#n
! 0 < x < a/2=1, 0 < y < b=1, Ne= 20x20 (dx=dy=0.05)
! Dxx = k^2*Dyy, k=2
! generated by MeshGenRect v54.3 of 22/08/18
--------------------------------------------------------------------------------
! parameters:
Ne	! total # of elements
Nn	! total # of nodes
Nb	! total # of BC faces
Nm	! total # of materials
Nd	! max elements per node
Ng	! max Gauss points in each coordinate

! properties:
Sx	! cross section area normal to x

! biCGstab input:
! ipar(1:3), ipar(5:6), fpar(1:2), fpar(11) - see comments in biCGstab
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! initial setting TEMPORARY, to be done in READIN or INIT
	ipar( 1) = 0	! start bCGstab solver
	ipar( 2) = 0	! no preconditioning
	ipar( 3) = 0	! default stopping criteria
	ipar( 4) = ldw	! work array dimension for bCGstab
	ipar( 5) = 0	! Krylov subspace dimension (unused for bCGstab)
	ipar( 6) =-1	! max. [A]{v} operations (unlimited if <0)
	ipar(11) = 0	! initialize FLOP counter
	fpar( 1) = tol	! relative tolerance
	fpar( 2) = tol	! absolute tolerance
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
--------------------------------------------------------------------------------
formers:
--------------------------------------------------------------------------------

